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Abstract 

We describe a unitary matrix model which is constructed from discrete analogs of the 
usual projective modules over the noncommutative torus and use it to construct a lat- 
tice version of noncommutative gauge theory. The model is a discretization of the non- 
commutative gauge theories that arise from toroidal compactification of Matrix theory 
and it includes a recent proposal for a non-perturbative definition of noncommutative 
Yang-Mills theory in terms of twisted reduced models. The model is interpreted as a 
manifestly star-gauge invariant lattice formulation of noncommutative gauge theory, 
which reduces to ordinary Wilson lattice gauge theory for particular choices of param- 
eters. It possesses a continuum limit which maintains both finite spacetime volume 
and finite noncommutativity scale. We show how the matrix model may be used for 
studying the properties of noncommutative gauge theory. 
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1 Introduction 



Noncommutative gauge theory first appeared in string theory within the framework of 
toroidal compactifications of Matrix theory [1]. It was argued that compactification on 
a noncommutative torus corresponds in 11-dimensional supergravity to null tori with a non- 
vanishing light-like component of the background three-form tensor field. Subsequently it 
was realized that the deformation to a noncommutative space can be described in Type II 
string theory as the effect of turning on a constant Neveu-Schwarz two-form tensor field B^ u 
in the worldvolumes of D-branes [2, 3]. The parameter 9 which deforms the space of func- 
tions on the worldvolume to a noncommutative algebra is related to the 5-field background 
by 9 ~ B~ x . The low-energy effective field theory for the gauge fields living on the D-brane 
worldvolume can be described by a noncommutative gauge theory. 

A non-trivial issue concerns the renormalizability of such gauge theories, given their 
unusual non-polynomial interactions. The perturbative renormalization properties of non- 
commutative Yang-Mills theory have been studied in [4]. In this paper we will present 
a constructive definition of noncommutative Yang-Mills theory which is the analog of the 
usual Wilson lattice gauge theory [5] in the commutative case. Such a model has the po- 
tential of clarifying issues of renormalization as well as shedding light on non-perturbative 
aspects of the gauge theory. A concrete definition of noncommutative gauge theory has been 
proposed recently in [6], and further studied in [7, 8], based on a large N reduced model 
[9]— [11]. In this case an ultraviolet regularization is naturally introduced at finite N and is 
removed in the large N limit with an appropriate fine-tuning of the gauge coupling constant. 
One expects the resulting theory in the continuum limit to have three scale parameters, the 
extent L of the space-time, the scale A of noncommutativity, and the usual gauge theoretic 
scale parameter A. However, it is found that j scales as \/~N, which means that one is 
inevitably led either to a finite L with A = (commutative finite space) or to a finite A with 
L = oo (noncommutative infinite space) [7]. In the following we will show that there exists 
a more general constructive definition of noncommutative gauge theory which possesses a 
continuum limit whereby noncommutativity is compatible with a finite volume space. 

The noncommutative gauge theory that naturally arises from toroidal compactification of 
Matrix Theory [1] comes from the matrix model which is obtained by dimensionally reducing 
ordinary Yang-Mills theory to a point [11]. The action is 



where X^, \x = 1, . . . , d, are N x N hermitian matrices and d is the dimension of spacetime. 
To describe the compactification of this model on, say, a two-torus of radii Ri and R 2 , one 
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needs to restrict the action (1.1) to those matrices that remain in the same gauge orbit 
after a shift by 1^ in the direction ji. This is tantamount to finding configurations 

for which there exists unitary matrices /i — 1, 2, which generate the quotient conditions 
[1, 12, 13] 

X M + 21x^5^ 1 N = n u X^ v . (1.2) 

Taking the trace of both sides of this condition shows that these equations cannot be solved 
by finite-dimensional matrices. It is, however, straightforward to solve them by self-adjoint 
operators on an infinite-dimensional Hilbert space. The basic observation [1] is that consis- 
tency of the conditions (1.2) when represented on this Hilbert space requires that 

Q 1 Q 2 = e 2 ™ 9 n 2 Qi (1.3) 

for some real number 6. This means that the operators fl^ generate the noncommutative 
two-torus with noncommutativity parameter 9. 

The operators fi M may be represented on a Hilbert space T~L p ,q of functions fk(s) where 
sGR and fcGZ, [1]. By introducing a fixed, fiducial derivation V M on this Hilbert space 
which satisfies 

[V^O,] = 2^,0, , (1.4) 
a generic solution of (1.2) may be taken to be the sum of and a fluctuating part, 

X ll = iR v 5 l JV v + A ll {Z) (1.5) 

where generate the algebra of operators which commute with the O^'s. There is a standard 
construction of these operators on the Hilbert space H. p , q [1]. In noncommutative geometry 
this simply corresponds to the algebraic construction of vector bundles over the noncommu- 
tative torus and the solutions (1.5) are just connections on these bundles [14]. The bundle 
T~C Pt q is characterized by its "commutative" rank p = dimTi, pq \g = o and its magnetic flux q 
which is taken to be the constant curvature of the fixed connection V^, 2irq = Tri [V^, V u }. 
The gauge fields A^(Z) are then functions on a dual noncommutative torus and the substitu- 
tion of (1.5) back into the action (1.1) gives Yang-Mills theory on this dual noncommutative 
torus. 

This construction has been reinterpreted recently in terms of open string quantization in 
the presence of a constant background B- field [3]. The modules T~t p , q are constructed from 
the boundary worldsheet theory appropriate to one end of an open string terminating on a 
D2-brane and the other end on a configuration of p coincident D2-branes carrying q units 
of DO-brane charge. In this paper we will present a construction which is a straightforward 
discretization of the above formalism in terms of an N x iV unitary matrix model. We shall 
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recover all the parameters labeling the continuum theory in the large iV limit. In particular, 
in this formulation finite A is compatible with finite L. The proposal in [6] can be regarded 
as a special case, from which it becomes transparent why finite A is not compatible with 
finite L in that instance. Thus the ensuing matrix model naturally interpolates between 
the model in [6] and the continuum formalism in [1] for the matrix model of M-theory 
[12]. We will show that our model can be interpreted as a manifestly star-gauge invariant 
lattice formulation of noncommutative gauge theory, which reduces to Wilson's lattice gauge 
theory [5] for particular choices of the parameters even at finite N. We shall also describe 
how various aspects of noncommutative gauge theory can be systematically studied within 
the matrix model formalism. 

2 The Unitary Matrix Model 

We will describe the construction in the simplest two-dimensional case, but the generalization 
to arbitrary even dimension is straightforward. The model we consider is just a twisted 
Eguchi-Kawai model [9, 10], but with a certain constraint imposed on the matrices. The 
action is 

S = -f3j2 Z ^ t * (UJJJJpt) (2.1) 

where (/i = 1, 2) are N x N unitary matrices and Z^ = Z*^ is a phase factor called the 
"twist" . The constraint we impose on the matrices is 

n v Uvfil = e 277i5 ^ /N (2.2) 

where r M (fj, = 1,2) are integers which we will specify below. The constraint (2.2) is the 
exponentiation of Eq. (1.2). Now, however, the only condition required is that the U^s are 
traceless unitary matrices. It therefore represents a finite dimensional version of the quotient 
conditions for Matrix theory. 

To solve the consistency conditions (1.3), we take the unitary matrices 

= (r 2 ) m ® (f^ , n 2 = (i\r ® (f 2 ) t , (2.3) 

where T M and T M are unitary matrices of dimension M and q, respectively, which satisfy the 
Weyl-'t Hooft commutation relations 

r\r 2 = e 2wi/M r 2 r 1 , f ± r 2 = e 27Ti/q f 2 f 1 . (2.4) 

These algebras can be represented by the shift and clock matrices (I^),^ = Sj+i^, (r 2 ) jfe = 
^ e 2m/My-i fi.^ anc j s i m ii ar iy f or the T^. The integers M and q satisfy iV = Mq and we take 



4 



M = mnq. The deformation parameter 6 of Eq. (1.3) is given by 

e = *-™. (2.5) 

q nq 

The incorporation of two independent integers m and n in the above will enable us to take a 
large- N limit whereby the appropriate continuum limit is reproduced. They play a certain 
dual role to one another as we shall see. We will also see later on that the meaning of the 
integers p and q will be the same as in the description of the modules Ti Piq above. 

Our first task is to solve the constraint (2.2). We take the simplest particular solution 
Ufj, = associated with 1 r 1 = r 2 = mq: 

D 1 = (r\)t ® l q , D 2 = Y 2 ®l q . (2.6) 

These operators will become fixed covariant derivatives in the continuum limit. We then 
decompose using as 

U, = %D, , (2.7) 
where are unitary matrices which satisfy the constraint 

n u %nl = % . (2.8) 

These constrained matrices will become the gauge fields of the model in the continuum limit 
and they can be constructed as follows. Assuming that p and q are co-prime, we choose 
integers a and b such that 

ap + bq = 1 . (2.9) 

We then introduce unitary matrices 

z 1 = (r 2 ) n ®(r 1 y , z 2 = (r^ n (g) (f 2 y , (2.10) 

which commute with Q^. The commutation relation of the is 

Z 1 Z 2 = e 2md ' Z 2 Z X , (2.11) 

where 6' is given by 

6' = ^- a - (2.12) 
mq q 

and it is related to 9 through the discrete Mobius transformation 

* = (213) 

p — qu 

1 One can consider a more general particular solution = (D /1 ) 1 ' 1 where 1^ are integers. If is a divisor 
of mq, then this solution will reproduce in the continuum limit the noncommutative gauge theory associated 
with a torus of modulus R1/R2 = h/h- 
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In the continuum, the transformation law (2.13) would be just that between Morita equiv- 
alent noncommutative tori [14]. In fact, by identifying 6 with a constant Neveu-Schwarz 
two-form field, it is just the T-duality transformation rule for the 5-field [3, 15]. The rela- 
tionship between Morita equivalence and duality [15] means that certain noncommutative 
gauge theories are physically equivalent to one another. We will return to this point later 
on. 

Using Z M , we can define a basis for the solution space of (2.8) as 

J mi ,m 2 = (Z 2 )^(Z 1 )^ e™ e ' mim2 (2.14) 
where the phase factor is included so that 

J—mi,—m2 (Jmi,m2^ ■ (2-15) 

Since (Z^) mq = ljv, J miim2 is periodic with respect to m\ and m 2 with period mq. We can 
therefore restrict the integers m\ and m 2 to run from to mq — 1. It will prove convenient to 
introduce a lattice with (mq) 2 sites on the torus and to instead work with the basis defined 
by 

A(x) = J ™ e~ 2m ^ m ^ IL , (2.16) 

mi,m 2 

where x^ — 0, e, . . . , e{mq — 1) belongs to the lattice of the spacing e and the extent of the 
lattice is 

L = emq. (2.17) 
We have defined A(x) in such a way that the identities 

D„A(x)Dl = A(x-efi) (2.19) 

hold. Here efi denotes a shift by e of x M only. The relation (2.19) expresses the identification 
of the matrices _D M as discrete covariant derivatives. Note also that A(x) is hermitian, due to 
(2.15), and is periodic with respect to x\ and x 2 with period L. The proof of completeness 
the generators (2.16) is given in Appendix A. Given this complete set of solutions, we can 
write a general solution to (2.8) as 

^ = ^-2zZ U ^ A ^- ( 2 - 2 °) 
Using orthogonality we can invert (2.20) to give (see Appendix A) 

U,(x) = ±tr(u,A(x)) . (2.21) 
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In order that the right-hand side of Eq. (2.20) is a unitary matrix, the coefficients U^x) 
should satisfy a certain condition which will be given below. 

Having solved the constraint, our next task is to rewrite our model entirely in terms of 
Un(x), which are the gauge fields of the theory. For this, we use the identity (2.21) to regard 
A(x) as a map from the space of N x N matrices which commute with to the space of 
fields on a periodic L x L lattice. We use the natural definition of a product of two lattice 
fields fi(x) and f2(x): 

h(x)*f2(x) = Itr^/aACr)) , (2.22) 

where fi are the N x N matrices defined by f\ = (mq)~ 2 S ^ JX fi(x)A(x). This product is 
associative but not commutative. One can write it explicitly in terms of fi(x) as 

fi(x) * f 2 (x) = — ^ J2 h(v)M z ) e 2iB ^-y^ x »-^ (2.23) 

y,z 

where 

B — ———^ . (2.24) 
9'L 2 V ; 

These formulas are similar to Ref. [8]. (See also [16] for earlier works in this regard.) The 
product (2.23) can be considered as the lattice version of the star product in noncommutative 
geometry. To see this, we note that in the continuum the star product of two functions fi(x) 
and f 2 (x) may be defined as 

fi(x)*f 2 (x) = h(x) exp(^V^) h{x) ■ (2.25) 
Using a Fourier transformation, this definition can be turned into an integral form 

fi{x)*f 2 (x) = JJdydzK(x-y,x-z)f 1 (y)f 2 (z) (2.26) 

where the integration kernel K is given by 

K(x -y,x-z) = -jr— e-^-^^-v^-z.) _ ^.27) 

7r | det Ofxv | 

The expression (2.23) can be obtained from (2.26) just by restricting the variables x, y, z 
to run over lattice points. In this sense, the product (2.23) is a natural lattice counterpart 
of the star product in the continuum. We shall therefore call (2.23) a star product in what 
follows. 

Using the star product, we can write down the condition on U^ix) which is required for 
Ufj, to be unitary as 

U^x)* *Up(x) = 1 . (2.28) 
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In other words, the lattice fields U^x) must be star-unitary. We may now rewrite the action 
(2.1) as 



x 



X [iyiv 

^ Up(x)*U v (x + eti)*U*(x + ei>)*lC(x) , (2.29) 

Tfl 

X fM^tU 



where Z^ = Z^ u e _27 ™W M can be considered as a background rank-two tensor field. One 
can make Z^ v = 1 by choosing Z^ v = e 27 ™W M . Then the vacuum configuration is given by 
= In, which corresponds to U^ix) = 1, up to the symmetry of the model which we now 
proceed to discuss. 

The action (2.1) and the Haar integration measure for the matrices are invariant 
under the SU(iV) transformations 

Up^gU^. (2.30) 

The constraint (2.2) in general breaks this symmetry down to a subgroup of SU(iV). However, 
the constrained model is still invariant under (2.30) for any g that commutes with Q^. We 
can represent such a g in terms of a function g(x) on the lattice as 

where g(x) should satisfy g(x)* * g(x) = 1, but is otherwise arbitrary. The transformation 
(2.30) can now be interpreted in terms of U^(x) as 

U^x) ^ g(x)*U^x)*g*(x + efi) . (2.32) 

Therefore, the resulting theory of the lattice field U^(x) is manifestly invariant under this 
star-gauge symmetry. 

We can show that the theory (2.29) reduces to Wilson's lattice gauge theory [5] for 
particular choices of the parameters. Note that we can always make 9' = by taking n = ma. 
In this case, the star product becomes the ordinary product of functions. Therefore, U^x) 
becomes an ordinary U(l) field on the lattice and the action (2.29) becomes the ordinary 
Wilson plaquette action. We can also show that the integration measure for U^(x) is actually 
the Haar measure for integration over the group U(l) 2(m ^ 2 . Note that the Haar measure 
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for the N x N matrices and the constraint (2.2) are invariant under — > gU^ for any 
g which commutes with This can be translated into the invariance of the integration 
measure for U^x) under U^x) — > g^U^x). The uniqueness of a measure with such an 
invariance proves our statement. Thus, our lattice formulation of noncommutative gauge 
theory includes Wilson's lattice gauge theory on a periodic lattice of finite extent as the 
9' = case, even at finite N. We remark that in this case, although the matrices can be 
taken to be diagonal, the (mq) 2 degrees of freedom of the lattice gauge theory are contained 
in the N = mnq 2 = a{mq) 2 diagonal elements of U^. 

Going back to the general case of arbitrary 9', let us now consider the continuum limit 
of the model (2.29) when the lattice spacing e — > 0. We introduce the continuum field 
and operator through 

% = e ieA » , = e ied » . (2.33) 

The large N limit dictated by the continuum theory [1] is m ~ n ~ \fN and e ~ l/y/N 
with fixed a, b, p and q. Both L given by (2.17) and B given by (2.24) are finite in such 
a large N limit. The resulting gauge theory is constructed from connections of a rank 
p bundle of magnetic flux q. We will see this explicitly in the next section. Note that 
the field theory (2.29) is actually of rank 1. This is one of the characteristic features of 
Morita equivalence or alternatively of T-duality transformations between different brane 
worldvolume field theories. The original SU(p) Yang-Mills theory on the noncommutative 
torus with deformation parameter 9 is physically equivalent to a U(l) Yang-Mills theory on a 
dual torus with noncommutativity parameter (2.13) that implicitly contains the information 
about the rank p of the underlying vector bundle. The case q = 0, representing a trivial 
gauge bundle, can also be constructed and will be presented elsewhere. 

However, as far as the continuum limit of the lattice theory is concerned, we need only 
send m to infinity, but not necessarily n. If n is finite as m — > oo, this does not lead 
to the solutions constructed in the continuum [1] for hermitian operators and is instead 
associated with unitary operators acting on periodic functions of < s < nq. The particular 
case of q — n — 1, for which the condition (2.2) is trivial and our model reduces to the 
ordinary, unconstrained twisted Eguchi-Kawai model, is of this type. It corresponds to the 
interpretation of the twisted large N reduced model in terms of noncommutative gauge 
theory which was proposed in [6]. Since 9' = 1/N in that case, in order to make B finite 
one needs e ~ l/y/~N, which inevitably makes the physical extent of the torus scale as 
L ~ VN, reproducing the observation made in [7]. Note that the issue of whether or not 
a continuum limit really exists is a dynamical question that can be addressed, for example, 
by Monte Carlo simulation. A numerical simulation of the two-dimensional Eguchi-Kawai 
model has been done in [17], where a non-trivial large iV scaling behavior was found with 
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the parameter Ne 2 fixed, which is exactly the large N limit required to make the physical 
scale of noncommutativity finite. This in itself means that noncommutative gauge theory 
with a background tensor field can be constructively defined. 

3 Observables of Noncommutative Gauge Theory 

We will now describe how the properties of noncommutative gauge theory can be completely 
reformulated in the language of the unitary matrix model above. Let us define a lattice 
path which consists of s links by C = {/ti, . . . , fi s } and = {fi s , . . . , jii} for an opposite 
orientation. The path C connects lattice sites separated by the vector ^ = while 

3=1 

belongs to C. We introduce the following products of matrices along the path 

s 

D{C) = IJD W , DiC- 1 ) = D(Cy ; 

3=1 

u{c) = n (^a) • 

3=1 

Given the property (2.19) we then have 

A(x + £)= D(C)A(x)D(C- 1 ) , (3.3) 

where the right-hand side is path-independent because of the properties of the _D M . This 
results in the formula 

1 tr (AA{x)) * 1 tr (B A(x + £)) = 1 tr (AD(C)5D(C'- 1 ) A{x)) , (3.4) 

provided that A and B belong to the commutant of the algebra generated by Using 
(3.4), we can construct the matrix analog of the noncommutative phase factor along the 
lattice path which defines parallel transport for the gauge bundle in the continuum limit, 

U(x; C) d = f * f[ U h ( x + fc-i) = ^ tr (U(C)D(C- 1 ) A(x)) (3.5) 

3=1 

where the product in the middle is the star product. Under the SU(iV) gauge transformation 
(2.30) where U(C) — > gU(C)g J( , the right-hand side of Eq. (3.5) transforms as 

U{x- C) -+ ^ tr (gU(C)g^D(C- 1 ) A(x)) = g{x) *U(x; C) *g*(x + i) (3.6) 

10 



(3.1) 



(3.2) 



as it should for the phase factor. This formula extends (2.32) to an arbitrary open path. 

The continuum limit of the above construction is given by the large- N limit of the matrix 
model. We introduce d£ M = e/t, so that Eq. (3.1) takes the form 

e = J de , (3.7) 

and write down the continuum analogs of Eqs. (3.2) and (3.5) using (2.33) as 

D{C) = P exp (i J de d„) , 

U{C) = P exp (i J de (4 + <g) (3.8) 

and 

U(x; C) = * J] (l + i de A,(x + 0) • (3.9) 

Here we have defined the field A(x) by 

U^x) =*e ieA »W (3.10) 

where the exponential is understood as a power series expansion with the star-product. 
Notice that the in Eq. (3.8) cannot be absorbed by a shift of since does not 
commute with the Q's. This is the difference between the present model and the continuum 
limit of the twisted Eguchi-Kawai model where this translation is usually done. 

The phase factors (3.9) can be used to define a new class of observables in the matrix 
model, associated with noncommutative gauge theory. The standard closed Wilson loops 
W(C) of twisted reduced models [10] which are invariant under (2.30) can be expressed via 
U(x; C) as 

W(C) = 1 tvDiC- 1 ) itr U(C) 

= ^E^rKWC-'lAW)^^^), (3.1!) 

since D(C~ V ) is a c-number. Therefore, the analog of W(C) in noncommutative gauge theory 
is a sum over lattice points of U (x; C) , which is understood as the sum over translations of the 
closed path that preserve its shape. This object is star-gauge invariant due to this additional 
summation. For the simplest closed loop, i.e. the plaquette, it is used in constructing the 
action (2.1). What is rather surprising in noncommutative gauge theory is that one can 
actually construct a star-gauge invariant observable associated with an open path, as has 
been found in [7] . We will now describe how such observables appear in our model and point 
out an interesting consequence of the finiteness of the spacetime extent. 
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Star-gauge invariant quantities can be constructed out of (3.5) with the aid of a lattice 
function Se(x) which has the property 

S e (x) * g(x) * Se(x)- 1 = g(x + 1) (3.12) 

for arbitrary star-unitary functions g(x). Here again £^ is the relative separation vector 
between the two ends of the open loop. Star-gauge invariant quantities can then be defined 
by (mq)~ 2 J2 X ^e(x) -kU(x; C). The property (3.12) in the matrix model becomes 

S e A(x)Sl 1 = A(x - I) (3.13) 

where 

s '^E^w a w ( 3 - 14 ) 

belongs to the commutant of the algebra generated by Using the definition (2.16) and 
with a little algebra, we obtain that Se(x) should satisfy 

S t (x + 6'Ljl) = e- 27Tie ^ /L S e (x) (3.15) 

where L is given by (2.17). Assuming that 9'mq = n — ma and mq are co-prime, the only 
solution to (3.15) is 

S e (x) = e 2mk - x/L (3.16) 

where = 0, 1, • • • , (mq — 1) and 

i^O'Le^K + n^L (3.17) 

with an integer vector n M . As is seen from (3.16), the ratio l-nk^jL plays the role of the 
momentum variable and it is related to the distance vector by Eq. (3.17). The longer the 
open loop is, the larger the momentum 27rk fJi / L should be. The discretization of momentum 
due to the finite extent of the torus leads us to an interesting consequence that should 
also be discrete. In the commutative case when 6' = 0, we obtain £^ = n^L reproducing the 
known fact that the only such gauge invariant quantities are the Polyakov loops (holonomies 
of noncontractable loops on the torus). It is remarkable that in noncommutative gauge 
theory on a finite volume there exist other objects of this kind with discretized values of 
the distance It remains discrete in the continuum limit since L is finite. This is the 
difference from the analogous quantities constructed in [7] for the IIB model, where can 
be an arbitrary vector in the large iV limit. The matrix description of the star-gauge invariant 
open loop is given by 

= ^(uMDiC-^St) (3.18) 
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where St = Jk 2 ,-k 1 is given by (2.14) for the solution (3.16). Its star-gauge invariance can 
be directly checked by noting that D(C~ 1 )Se in (3.18) belongs for £ given by Eq. (3.17) to 
the commutant of the algebra generated by Z M , i.e. commutes with g. 

The matrix model determines the dynamics of noncommutative gauge theory. Let us 
demonstrate how the classical equation of motion emerges in the matrix language. For 
simplicity we take the continuum limit using the relation (2.33). The continuum action 
reads 

S[A] = tr ( {P^ - ) + tr (a^ ) (3.19) 

where 

Ffiv = id^A u — id v Ap, + % A^A V . (3.20) 
Here f^ u is the constant curvature of the gauge bundle given by 

-% [d^ d v ] = fp U ljv (3.21) 

where in the two dimensional case 

fuv = < ^—p i R\R2e llv . (3.22) 

p — qO 

In the construction of section 2, R x = R 2 = l/enq are the radii of the two-torus (see 
footnote 1). Eq. (3.22) is the standard formula for the curvature of the module H Ptq [1]. 
It should be understood, however, as being multiplied by the identity operator with 
Tr l PjQ = p — q9, so that the integral curvature of the bundle is Tr fp l/ /(2irRiR 2 ) = q. 
In the present case this trace operation corresponds to multiplying the curvature (3.22) 
by the dimensionless area factor a/ (.Rii^X-^i-^) = m/n giving the volume of a unit cell 
in the "phase space" of the d^s. This is analogous to the derivation of the dimension 
of the Hilbert space H P:q presented in [1]. The (infinite) hermitian matrices in (3.19) are 
unconstrained while the constraints are taken into account by the Lagrange multipliers a^ u . 
The action (3.19) is of the type considered in [18], but now with the additional constraints 
imposed on A. 

The variational derivative 

5 



■A v (y)=5 lu ,8(x-y) (3.23) 



can be represented in the matrix language as follows. Given (2.20), (2.33) and (3.10), we 
have 

A li = J dx A(x)A(x) (3.24) 

which implies 

tr(^^A{x)) . (3.25) 



SAp(x) \dAp 
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Equation (3.23) is now reproduced as 



tr 



d 



OA,, 



-A(x) ) -ti(A v A(y)) = ^ - tr ( A(x)A(y) ) = 5 fiu 5(x-y) 



N 



(3.26) 



as it should. We can treat the matrix elements of A u as independent because of the com- 
pleteness of the generators of the commutant. Notice that there is an ordinary product of 
the traces in (3.26) rather than the star product. Acting by (3.25) on the action (3.19), we 
get 

5 



SAJx 



S = tr 



d M + A„, - A(x) + tr cv [A(x), ^ 



tr ( [c^ + i M , F^ A (a;)) 



(3.27) 



Eq. (3.27) reproduces the noncommutative Maxwell equation. 

The matrix representation (3.25) of the variational derivative is actually most useful for 
deriving the Schwinger-Dyson equations of the quantum noncommutative theory and, in 
particular, the loop equations. To illustrate the technique, let us first calculate how the 
variation 5/5A^{z) acts on the noncommutative phase factor U(x; C), which determines the 
contact term in the loop equation [19]. Using (3.25), we get 



5 



6A„(z) 



U(x;C) = tr 



d 



OA. 



A{Z) ) ^(uWDiC-^Aixj) 



ft i 

= i J dC j^tr (uiC,) A(z)U(C 2 ) D(C- l )A(x)) 

= i [ d?U(x;C 1 )*6(x + Z-z)*U(x + Z;C 2 ) 
Jo 



(3.28) 



where C 1 and C 2 are the parts of the contour C, C = CiC 2 , separated by £. We can similarly 
calculate how the area operator d IM 5/5a^ u (z) (z G C) acts on U(x;C). This calculation is 
purely geometrical and gives 

d, U(x; C) = i tr {u{C x ) [d„ + A„ F, u ] [/(^(C*- 1 ) A(x)) 

= -iU{x-C l )^{d^^ + iA^J r ^-iJ r ^^A^{z)^U{z-C 2 ) (3.29) 



where 



!F d^Ay dyA^ ~\~ iA.^i ^ A. v iAy ^ A^. (3.30) 



That is, the operator d^S /Sa' J ' u (z) inserts the Maxwell equation in the noncommutative phase 
factor at the point z, as anticipated. 
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The standard loop equation of large- N Yang-Mills theory for the Wilson loop average 
(W(C)) emerges from Eqs. (3.28) and (3.29) in the q — n — 1 case by putting z — x, 
summing over all x and using the formula 

^2 ^ ij (x)A k \x) =N5 u 5 kj (q = n = l case) . (3.31) 

X 

This equation is quadratic in (W(C)) due to large-A^ factorization of correlators. To better 
understand the consequences of Eqs. (3.28) and (3.29) for q,n ^ 1, let us consider the case 
of 6' = 0, whereby the continuum limit is rank one commutative gauge theory, previously 
known as Maxwell's theory. The phase factor (3.9) for a closed loop is now gauge invariant 
since the star-product becomes the ordinary product, so that U(x;C) becomes the usual 
phase factor of electrodynamics which is independent of x while the g's cancel on the right- 
hand side of Eq. (3.6). The loop equation for the average of the phase factor can be obtained 
by combining the averages of Eqs. (3.28) and (3.29) which results in the standard linear loop 
equation 

s * jJw> ( u(c) ) = w I ie 5{i - z) ( u{c) ) ■ (3 ' 32) 

We have just illustrated by this simple example how the phase factors (3.9) can indeed 
correspond to observables in noncommutative gauge theories associated with the unitary 
matrix model. This is precisely the novel feature of the present matrix model that was 
pointed out in section 2, namely that in the large N limit it is possible to arrive at a U(l) 
continuum gauge theory. 
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Appendix A Proof of Completeness 

We will demonstrate that the A(x) defined by Eq. (2.16) form a complete set for the space 
of solutions to the constraints (2.8), i.e. that any N x N complex matrix A that commutes 
with f2 M (/i = 1,2) can be written uniquely as 

A = -±- 2 TA(x)A(x) . (A.l) 
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We first note that 

£ = [Aeg\{N,C) | An il = %A , /i = l,2} (A.2) 

defines a linear subspace of gl(iV,C) which has an inner product defined by tr (A^B) for 
A, B e gl(iV, C). The A(x) e £ form an orthogonal set, 

ltr(A(x)A(y)) = (mg) 2 ^. (A.3) 

We now consider the linear subspace £' of £ spanned by A(x). We wish to show that £' = £. 
For this, we introduce a convenient orthogonal basis of gl(iV,C). Define A(x) by 

A(x) = W" 1 W" 2 e - wiemim2 e 27r ^ m ^" , (A.4) 

mi,m.2 

where x runs from to nq — 1 and we put e = 1 for simplicity. These matrices commute 
with Zn, they are mutually orthogonal, and they satisfy Q^A^Q^ = A(x — /}). We take 
A(x)A(x) as an orthogonal basis of gl(AT,C). We now consider a generic element which 
belongs to the orthogonal complement of £' in gl(AT, C) given by 

£ f(x,x)A(x)A(x) . (A.5) 

Requiring that it commutes with both Qi and f2 2 implies immediately that f(x,x) = 0, 
which completes the proof. Using the orthogonality (A.3) of the basis A(x), we can write 
the A(x) in (A.l) as 

A(x) = ^tr jAA(i)) . (A.6) 
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